Non-Fermi liquid behavior and superconductivity in the dissipative model 
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We study an exactly-solvable model which shows a zero-temperature transition from a non- Fermi 
liquid to a Fermi liquid as a function of particle density. The quantum critical point separating 
these two states is not associated with the usual ordering transition of a bosonic degree of freedom. 
We analyze the properties of this critical point in terms of an effective density of states, and we use 
this concept to generalize this transition to a class of related quantum critical points. We study the 
superconducting instability of the low-density (non-Fermi liquid) regime and we show that it has an 
enhanced pair susceptibility relative to the Fermi liquid. Wc comment on the applicability of our 
results to low-carrier density superconductors. 



In recent years there has been great interest in un- 
derstanding the physics of metallic states which are not 
Fermi liquids — motivated largely by a host of experi- 
mental observations Q. One route to non- Fermi liq- 
uid behavior stems from proximity to a quantum criti- 
cal point 1^, ^ , in other words a zero-temperature phase 
transition which occurs as a parameter such as pres- 
sure or chemical composition is varied. Examples include 
CePd2Si2 [D , and CeCue-xAua; § which show quantum 
critical behavior as pressure and gold concentration are 
varied respectively. In both of these materials the criti- 
cal point separates a magnetically ordered state from a 
paramagnetic one. 

The metallic state of the cuprates also shows non- 
Fermi liquid properties. Furthermore there is evidence 
of a pseudogap energy scale , T* , which is suppressed 
to zero temperature at some critical doping pcrit , and 
it has been suggested Q that this forms a quantum crit- 
ical point albeit masked by superconductivity. However, 
if this quantum critical point is responsible for the non- 
Fermi liquid physics of the metallic state, it differs from 
the above examples in that it does not appear to have 
an associated ordered phase: thermodynamic measure- 
ments 1^ show that T* corresponds to a crossover and 
not a phase transition. It has been argued jl^ that there 
may be a hidden order parameter, but an alternative sce- 
nario is that the cuprates show quantum criticality with- 
out associated order at finite temperature. Such behavior 
is common in one dimension where thermal fluctuations 
destroy finite-temperature order, but this motivates the 
search for examples of zero-temperature phase transitions 
in d > 1 without any associated order. 

In this Letter we consider a model of dissipative 
fermions, first studied by Wheatley , which possesses 
a zero-temperature transition but with no order other 
than the Fermi surface itself. At low carrier density, this 
model is a non-Fermi liquid metal but with a zero tem- 
perature transition to a Fermi liquid at a critical density. 
We characterize the model's behavior at the critical point 
by the Sommerfeld coefficient and the paramagnetic sus- 
ceptibility as T — > 0. We present an interpretation of 
the phase diagram by introducing an effective density of 



states. This allows us to consider more general metal 
to metal transitions of this type. We also consider the 
superconducting instability of this model. Although our 
starting model is artificial, we speculate on the connec- 
tion between models of this type and atom traps contain- 
ing neutral fermions as well as low carrier concentration 
superconductors. 

Our starting point is to consider a single particle of 
mass M coupled to a bath of harmonic oscillators 
This bath has a spectral density J{uj) and was origi- 
nally introduced to mimic the dissipative effects of a 
real- world environment |p^ . However, our philosophy 
here is to introduce the oscillator bath to give the sin- 
gle particle a non-trivial (i.e. non-local) dynamics in 
imaginary time while retaining the exact solvability of 
the system. This may be seen on integrating out the os- 
cillators, which leads to a single-particle effective action 
Sc« = 5'ke + S'diss, where S'ke = \hM (3Y,n^^'^-^^^n is 
the kinetic energy and 
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Here r„ is the Fourier component of the path in imagi- 
nary time and Un = is a Matsubara frequency. One 
usually chooses J{uj) = rjuj with a high-frequency cut-off 
^Wc, SO that S'diss corresponds classically to a frictional 
force proportional to velocity, F — rjv. In the first part of 
this paper we also make this choice, but later we consider 
the consequences of more general spectral functions. 

An exact expression can be found for the single-particle 
partition function z{P) When the temperature is 

low compared with the damping rate, — 2M/ri, so 
TiP/tq :s> 1, this may be simplified to 
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in d-dimensions, where no ~ [M /2'K'hTi:,f^'^ has the di- 
mensions of density and Eq depends on the details of 
the cut-off in J{oj) and will be absorbed as a shift in 
the band-edge. The dissipation defines an energy scale 

Ed = Gh/dTTTQ. 
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Having defined the properties of the single particle we 
now consider a gas of such particles {i.e. fermions cou- 
pled to their own independent oscillator baths). Anti- 
symmetrizing appropriately for fermions is most easily 
done using a path integral approach which gives the 
grand thermodynamic potential for the dissipative gas 
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Wheatley Q solved this equation numerically to de- 
termine the chemical potential /i and the susceptibility 
X, finding a T = phase transition from a low density 
{n < no) phase with Curie susceptibility to a high den- 
sity (n > uq) phase with Pauli susceptibility. Here we 
consider the system properties in the vicinity of the crit- 
ical point {n no) and we find a generalization of the 
model. 

The chemical potential is given by the solution of 
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once we have determined /x, the formulas for x ^-nd 7 
follow straightforwardly. The solution of Eq. (jj) in the 
limit T depends on whether we approach the critical 
density from higher or lower densities. 

Approaching the critical density from above (d„ — 
1 > 0) we find 
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We see that proximity to the critical density sets a new 
scale EdSn- For low temperatures, PsdSn ^ 1, the prop- 
erties recover those of the free Fermi gas: 7 and x 
constant. We call this the Fermi liquid regime, although 
there are no Landau parameters since the particles do 
not interact with each other. 

Approaching the critical density from below, 5„ < 0, 
we find /i — kBThi{—S~^ — 1) and hence for small Sn 
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There is also a residual ground state entropy in the low 
density phase which is given by iS'(T = 0) = fcBno(|(5„| — 
\Sn\\n\S„\). 

It is rather straightforward to see how this new phase 
of matter comes about. Usually the transition to a quan- 
tum fluid occurs when the thermal wavepacket of a parti- 
cle starts to overlap with the wavepacket of neighboring 
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FIG. 1; (a) 7(r) for a range of densities, (b) T — n phase 
diagram of the dissipative gas. Solid lines are phase transi- 
tions, dashed lines are crossovers, SC = superconductor and 
the circle is a quantum critical point. 



particles. In this model dissipation cuts off the de Broglie 
wavelength such that for particle densities lower than no 
the quantum liquid regime is not entered at any temper- 
atures. For the case of a dissipative Bose gas one has 
the possibility of a Bose metal phase — the bosonic 
analogue of the non-Fermi liquid state that results here 
for dissipative fermions. 

Having characterized the properties in the vicinity of 
the critical density, we show in Fig. |l|a the results of a nu- 
merical calculation of the Sommerfeld coefficient as the 
critical density is approached. This clearly shows the di- 
vergence in 7. It also shows how the divergence is cut off 
at higher densities by the crossover into the Fermi liquid 
regime. We have calculated the form of the divergence 
at the critical point itself and find that x oc InT and 
7 oc In^T. At a given T the properties of the high-density 
regime below a crossover scale (5„ ~ T/ed qualitatively 
resemble those of the low-density regime, which we illus- 
trate by a dashed line on the phase diagram of Fig. ^3 
|l5| . There is a second, high-temperature crossover which 
is set by the highest scale in the oscillator bath /3ujc ~ 1, 
above which no dissipation is felt by the particles. 

Further insight into the nature of the phase diagram 
and the route to generalizing this model can be found 
by considering the "effective" density of states. We use 
this term advisedly for the energy eigenstates of a dissi- 
pative particle are not well defined — it continuously in- 
teracts with the oscillator bath. However, since we have 
a single-particle partition function, Eq. (||), we can use it 
to define an effective single-particle density of states ges 
via zil3) = dEgcs{E)e-'^'^. Thus the effective den- 
sity of states is formally the inverse Laplace transform 
of the single-particle partition function. From here it is 
easy to show that the thermodynamics can be obtained 
directly from gcs by direct substitution into the usual 
definition of the non-interacting thermodynamic poten- 
tial n ^ -T dE g{E)\n[l + e'^'^^'-'^^. The effective 
density of states therefore captures the thermodynam- 
ics of the dissipative gas, although it does not contain 
information about correlations. 
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For the partition function of Eq. (g) we find an effective 
density of states near the band-edge 
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It is the initial delta function in the effective density of 
states that leads to the non-Fermi liquid to Fermi liq- 
uid transition. The delta function can accommodate hq 
fermions so for n < rig the Fermi energy is pinned there 
and the entropy is clearly finite at T = 0; this is the 
non- Fermi liquid regime. When n > hq the Fermi energy 
moves a distance (5„ed above the delta function, to the 
region where geS is weakly energy-dependent, and the 
system behaves like a degenerate Fermi gas. However, 
when T ~ dn£d the Fermi function is sufficiently broad 
that the occupancy of the delta function fluctuates and 
there is a crossover to dissipative behavior. 

We now consider extensions of this approach to con- 
sider a wider range of Fermi liquid to non-Fermi liquid 
transitions that can be found for dissipative fermions. 
The obvious extension is to go beyond Ohmic dissipation 
and consider more general oscillator baths with power- 
law spectral densities J(w) — r/Lu°' up to a high-frequency 
cut-off ~Wc, where < a < 2. Integrating out the oscil- 
lators gives (say for d = 1) the partition function 
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At low temperatures, ksT <^ Utq ^ , this is approximately 
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The band edge shift Sq is again cut-off dependent. In 
terms of the energy scale Eq. = 6h/{2 — a)TTTo the inverse 
Laplace transform gives, to leading order in E/sa, 
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for a < 1 . 



So, for a < 1 an infinite number of particles can be ac- 
commodated in the ground state and the gas is a non- 
Fermi liquid for all particle densities. For a > 1 we 
always have conventional Fermi gas behavior. Thus only 
Ohmic dissipation gives rise to a quantum phase transi- 
tion between a Fermi liquid and a non-Fermi liquid. 

Nevertheless, the notion of an effective density of 
states give us an alternative route to explore other 
models. Since a quantum phase transition will al- 
ways be present when there is a delta function in g^s, 



we will work directly with a generalized gcs{E) = 
"o/ed [5{E/ed) + {E/e^Y 6(i?/£d)], rather than explic- 
itly choose the spectral function of the oscillator bath. 
We must have v > —1 for goff to be integrable. This 
gives the partition function 
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Substituting this into Eq. ^ and using the properties of 
the polylog Li^(a;) = X^n ^ we find at n = no 
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Thus only for a < is the quantum critical point char- 
acterized by divergences in x 7. Note that for all a 
the Wilson ratio 7/x is proportional to InT. 

We now address the applicability of this model to real 
physical systems. Our method of including all non-trivial 
effects via a bath of oscillators has parallels with dynam- 
ical mean field theory where interactions are included via 
a self-consistent frequency-dependent self-energy |l^, ^ . 
We have not sought a self-consistent approach since our 
goal was to categorize all the possible transitions in the 
model. Nevertheless our starting model may be directly 
relevant to physical systems such as low density Fermi 
gases which couple primarily to low energy bosonic modes 
rather than each other. In particular, this necessitates 
short range interactions. Examples could include neutral 
^Li ions in atom traps where the bosonic modes would be 
the "optical molasses" of laser cooling giving a dissipative 
environment JTot . 

In a condensed matter context we apply this to 
very low carrier density metals with strong electron- 
phonon coupling but in a temperature window where the 
Coulomb repulsion does not dominate (above any ten- 
dency to Wigner crystallize for example). One possible 
case is SrTiOa. This is a band insulator |20| but be- 
comes a superconductor with Tc < 1 K |^ on electron 
doping to a concentration of about 10^^ cm~^. The mech- 
anism of superconductivity in this material has long been 
in question not least because superconductivity only ap- 
pears below a critical density. It is therefore natural to 
ask whether the non-Fermi liquid state we have consid- 
ered is unstable to pairing. 

We calculate the pair susceptibility in the dissipative 
model via the pair propagator ||2^ . In d dimensions the 
real-space density matrix p{r, (3) for a dissipative particle 
is given approximately by 
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where ^ = •\/7r?iro/mln2 is a dissipation length. The 
Matsubara Green's function can then be calculated using 
G^^ — p{r, x) Qip-+^^^r,)x which follows from the 
definitions of G^)^ and p [|2^. This integral is only well 
defined for fJ. < 0, which is the case for densities n < 
but it can be analytically continued to all /i, as can be 
checked by evaluating it for free fermions. Substituting 
Eq. (Ill) with d=3 yields 



3/2 



2'!Th 



2n 1/2 



3/2 



1 



fi + ihuJn 



(17) 



From the definition of the pair propagator K{r) = 
T J2uj Gui„ (r) G-uj,^ (r) we can determine Tc in the usual 
way as the solution of 1 = J dr V{r)K{r, (i). We consider 
a BCS-like attractive potential with small finite width tq. 
For n < uq this leads to 
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where the normalized potential Vgff depends on the de- 
tails of V{r) at r < Tq. For n > no, on the other hand, 
Tc is heavily suppressed and superconductivity does not 
occur for arbitrarily small Vcg. The dependence of Tc on 
n is shown in Fig. |l]b. The pairing enhancement can be 
traced back to the delta function in gcs{E). 

However, we can expect a different field dependence of 
the upper critical field, Bc2(T), for superconductivity in 
the dissipative model because of the appearance of the 
dissipation length scale. A magnetic field destroys super- 
conductivity by dephasing the two electrons that make 
up a Cooper pair; in a dissipative gas the electrons are 
already dephased due to the oscillator bath and so we ex- 
pect a field to have negligible effect at small T. Bc2{T) is 
given by the solution of 1 = /^^^ dr K{r^ (3) e'^^"'^^ 1"^^ 
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no we find just below Tc that Bc2 = 
37T cTo \" T y ^^'^ temperatures Bc2{T) ~ InT. 

In practice this divergence in Bc2 [T) would be cut off by 
Landau quantization effects which we have not included. 

So we expect a significant enhancement in the up- 
per critical field in the dissipative model compared with 
Bc2(0)/T2 which for a Fermi liquid is ^n^^i^^Y- In 
thin films of SrTiOa this ratio is enhanced by about 
an order of magnitude relative to the Fermi liquid ex- 
pression, where m* = 5.3mo is determined from 7 [ p5| . 
Thus SrTiOa qualitatively resembles the dissipative gas 
both in its phase diagram and its relative insensitivity to 
applied magnetic fields. 

To summarize, we have considered the Fermi liquid to 
non- Fermi liquid zero-temperature transition in a model 
of dissipative fermions. This critical point occurs without 
any associated ordering and may be generalized via the 



concept of an effective density of states. At low temper- 
atures the non-Fermi liquid shows an enhanced pairing 
susceptibility and a phase diagram similar to the low den- 
sity superconductor SrTiOa . More generally, this type of 
non-Fermi liquid state could occur during laser cooling 
in fermionic atom traps. 
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